SPREADING FACTOR ESTIMATION SYSTEM AND METHOD 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application claims priority from provisional applications: Serial No. 
60/263,669, filed 01/23/01 . The following patent applications disclose related subject 

matter: Serial Nos. 09/ , filed .... ( — ). These referenced applications have a 

common assignee with the present application. 

BACKGROUND OF THE INVENTION 

The invention relates to electronic communications, and more particularly to 
CDMA-based coding, transmission, and decoding/synthesis methods and circuitry. 

Code division multiple access (CDMA) coding has been extensively used in 
such applications as cellular and satellite communications. CDMA signals increase 
f| the spectrum required for the transmission at a particular data rate by modulating 

each data symbol with a spreading code having a rate larger than the data rate. The 
H same spreading code is used for each data symbol. Typically, the spreading code 
y comprises of a few tens or a few hundreds of elements, called chips. To decrease 
J^. the correlations among spreading codes assigned to different users, and thereby 
nj reduce the Interference among different users, the data stream after spreading is 
J typically scrambled with a pseudo-noise (PN) code that is generated serially and 
3 cyclically and has a larger period than the spreading code. Examples of such CDMA 
signal spreading are the schemes used by the IS-95/CDMA2000 and 3GPP systems. 

With CDMA, the signals from all users simultaneously occupy the same 
frequency band, and the receiver separates the multiple signals by exploiting the 
crosscon-elation properties of the spreading and scrambling codes that are applied to 
each user's signal. The receiver attempts to match In time the spreading and 
scrambling codes of the desired signal with a replica of these codes. Only then the 
demodulation result is meaningful; otherwise it appears noise-like. Thus, if the 
arriving signals have different codes or different code offsets, they can be 
discriminated at the receiver. The CDMA code for each user is typically produced as 
the modulo-2 addition of a Walsh code with a pseudo-random code (two pseudo- 
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random codes for QPSK modulation) to improve the noise-like nature of the resulting 
signal. A cellular system cell as illustrated in Figure 1 could employ IS-95 or 3GPP 
for the air interface between the base station and the mobile user station. 

The 3GPP system employs a 3.84 MHz bandwidth (a chip duration of 260 ns), 
and the spreading code length applied to each data symbol may vary from 4 chips to 
256 chips. The number of chips per symbol is called the spreading factor (SF), and 
the SF stays constant during each 1 0 ms duration frame. However, the SF may 
change from frame to frame. Each frame is partitioned into 15 time slots (0.667 ms 
duration) with each time slot consisting of 2560 chips. Thus the number of data 
symbols per slot ranges from 10 to 640, and the number of data symbols per frame 
ranges from 150 to 9600. 

In the 3GPP frequency division duplex (FDD) mode uplink, the Dedicated 
Physical Control Channel (DPCCH) is always spread by the code Cc=Cch,256,o with SF 
= 256 while the Dedicated Physical Data Channel (DPDCH) is spread by the code 
Cd,i = Cch.sF.k where k= SF/4 is the code number. Note that 3GPP uplink uses QPSK 
modulation with single-channel DPDCH and DPCCH being in phase and quadrature, 
respectively. The recursion defining even-numbered channel codes is concatenation: 

Cch.2"*\2k = Cch,2" k . Cch/. k with initial code Cch.4,1 = "+1, +1,-1,-1". Thus a 
consequence of the code selection k = SF/4 is that for two codes with SF's of length 
SFi and SF2, where SFi < SF2, the part of the code for SF2 consisting of the first SFi 
chips is identical to the code for SFi. In the uplink of 3GPP, the SF may take the 
values 2', 1=2,3,4,5,6,7,8. Because the code used for different SF's is a repetition of 
the code for SF=4, despreading for all possible SF's can be combined in a single 
process and as a result the number of despreaders need not increase. 

The data rate and the SF used in each frame are specified by the Transport 
Format Control Indicator (TFCI) field that is transmitted in the DPCCH control 
channel. The DPCCH control channel consists of 15 time slots per frame, and the 
TFCI bits are time-multiplexed, together with other control information, among the 15 
slots. The DPCCH control channel is always spread by the all-zero code 

("+1 ,+1 ,+1 +1 ,+1") with SF=256 and can therefore be always despread once the 

received signal paths have been identified through a searching process at the 
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receiver. By despreading the DPCCH control channel, at the end of a frame the 
receiver knows all of the TFCI symbols and can determine the data rate and the SF 
of the DPDCH data channel. Therefore, in order to despread the DPDCH data 
channel, the received data need to be stored for one entire frame period before 
despreading. Once the SF becomes known, the stored data need to be quickly 
transferred to the baseband part of the receiver to reduce the associated processing 
delay. 

Thus the known method needs to store an entire frame prior to despreading 
and baseband processing; this requires large front-end memory, large data rate 
transfers between front-end memory and hardware and creates processing delay. 

SUMMARY OF THE INVENTION 

The present invention provides estimation of the spreading factor in less than 
p a frame by despreading with a plurality of spreading factors and estimating the actual 
5 spreading factor from the despreading results. 

f;f; The preferred embodiment SF estimation methods have advantages including 

y I 

yj substantially reducing the requirements of large memory, high data rate transfers, 
y and processing delay. Also, the preferred embodiment implementations of the 

HJ estimation methods only require functions (e.g., despreading and Rake combining) 

Q 

^ that already exist in an uplink receiver. 

S 
N 

BRIEF DESCRIPTION OF THE DRAWINGS 
The drawings are heuristic for clarity. 

Figure 1 is a diagram of a typical communications cell comprising of several 
users and a base station. 

Figure 2 is a block diagram of a conventional spread spectrum CDMA Rake 
receiver with despreading performed for every possible spreading factor. 

Figure 3 is a block diagram for the initial decision statistics after Rake 
combining of the preferred embodiment SF estimation methods. 

Figures 4.1-4.7 show a set of decision statistics for a firet prefen^ed 
embodiment SF estimation method. 
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Figures 5.1-5.7 show another set of decision statistics for a first preferred 
embodiment SF estimation method. 

Figures 6.1-6.7 show a set of decision statistics for a second preferred 
embodiment SF estimation method. 

Figures 7.1-7.7 show another set of decision statistics for a second preferred 
embodiment SF estimation method. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



1 . Overview 

Preferred embodiment methods estimate received CDMA signals' spreading 
factor (SF) from a plurality of despreadings with hypothesized SF's and decide on an 
SF after a specific observation period ensuring a sufficiently small error probability. 
Figures 2-3 show functional blocks for some preferred embodiment systems using 
preferred embodiment methods. The preferred embodiment methods may be 
implemented on programmable processors (such as DSP's) or in specialized circuits 
(ASIC's) or combinations such as system-on-a-chip with attendant analog-digital 
signal components. 

2. 3GPP preferred embodiments setup 

The first preferred embodiment methods and systems relate to the base 
station receiver and the uplink in a 3GPP code division multiple access (CDMA) 
communication system. Figure 1 illustrates a typical cell in a CDMA communication 
system where several mobile users communicate with a base station. 

For 3GPP the 10 ms frame period is divided into 15 time slots with 2560 chips 
per slot. There are ten Dedicated Physical Control Channel (DPCCH) symbols in 
each slot, and each symbol extends over 256 chips (spreading factor of 256). Those 
symbols may be for pilot, TFCI, power control, and feedback information for transmit 
antenna diversity. The DPCCH always has a spreading factor (SF) equal to 256 and 
uses the all-zero Walsh code. 

In contrast, the number of Dedicated Physical Data Channel (DPDCH) 
symbols in each slot depends upon the SF and may vary from 10 for SF=256 to 640 
for SF=4; although the spreading factor remains constant throughout a frame. That 
is, the 2560 chips of a slot may be partitioned into from 10 to 640 symbols with 256 to 
4 chips per symbol and yield corresponding data rates of from 15 to 960 ksymbols/s. 

Figure 2 schematically shows a first part of prefenred embodiment receiver 
circuitry used for SF estimation. For a particular user, let the total number of possible 
spreading factors (SF's) be denoted M and the number of received separable 
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multipaths be denoted L. Typically, M = 7 and L = 1-6. For each mobile user of 
unknown SF, the SF estimation circuitry consists of M Rake receivers, one for each 
of the possible M SF's. Only one of the M Rake receivers would exist if the SF were 
known for a user (and no SF estimation would be necessary for this user). The 
receiver includes a plurality of despreaders, a plurality of timers, a plurality of spread 
spectrum processors, and a plurality of adders. Users' signals are received at the 
antenna and subsequently downconverted and digitized. Since the various user SF's 
are not known, the received signal samples need to be stored in buffer memory until 
the SF is determined and despreading can occur. The searcher identifies the 
multipaths for each user. Subsequently, L identified paths for each user's signal are 
passed to the despreaders (groups of L fingers for each SF). Each despreader (or 
finger) uses a replica for a possible spreading code which is appropriately delayed in 
time to synchronize with the corresponding multipath. The 1^' despreader output for 
the 1®* path of the m*^ possible SF is denoted as Dm,i. Each despreader also includes 
a chip-code signal generator, a modulo-2 adder, and a summer. In addition to the 
despread signal paths, each despreader (possibly together with another unit) also 
provides a channel medium estimate for the corresponding path of a user's signal 
(Cm,i, Cm.L for the m*^ SF). The channel medium estimates are needed by the 
conventional Rake receiver to combine the L multipaths. In the 3GPP uplink, the 
channel medium estimate and the search may be based on the pilot symbols that are 
embedded in the DPCCH control channel; and each user having a differing 
scrambling code allows for the user separation. After despreading, the multipaths of 
the signal for each user and the corresponding channel medium estimates are 
provided to a Rake receiver for decision. The m ^ Rake receiver performs maximal 
ratio combining (MRC) on the multipath components of the signal and provides soft 
output Sm for the m * SF. It should be noted that the actual despreading process for 
all M possible SF's can be performed by a single despreading unit which can produce 
intermediate despreading results for the smaller SF's before completing the 
despreading process for the largest SF. This means that the total number of 
despreaders needed for all M possible SF's is equal to the number of multipaths L for 
the particular user and there is no increase in the hardware complexity associated 
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with the despreading operation. For clarity purposes, M replicas of such despreaders 
are shown in Figure 2. 

For the m"^ spreading factor, sum the absolute values of the Rake soft output 
Sm over a number of symbol periods (or over a number of slots); see Figure 3. The 
number of slots is selected so that the error probability for the SF estimation methods 
is substantially smaller than the frame error rate (FER) and comparable to the TFCI 
decoding error rate. Also, the number of slots required to achieve the desired 
performance for the SF estimation methods will determine the size of the memory 
required to store the received signal samples (Figure 2, buffer memory) and the 
ensuing data transfer rate and processing delay. Thus, the number of slots should 
be minimized, subject to the condition of achieving the desired performance for SF 
estimation, and in 3GPP it should not exceed 15 slots which is the frame period. 
After summation, the results may be divided by the square-root estimate of the noise 
and interference variance during the summation period; see the multipliers in Figure 
3. This is an optional operation that makes the values of the decision statistics 
dependent on the signal-to-noise ratio (SNR) as it will be later shown. Finally, the 
decision statistics are provided to the SF estimation (decision) unit for further 
processing and for a decision regarding the SF. This processing in the SF estimation 
unit defines the SF estimation method. Two preferred embodiment methods and the 
decision process are described in the following sections. 

In the 3GPP uplink, the channel medium estimate and the search may be 
based on the pilot symbols that are embedded in the Dedicated Physical Control 
Channel (DPCCH). The DPCCH always has a SF of 256 and uses the all-zero 
Walsh code. For that reason, the channel medium estimate will be the same for all of 
the M Rake receivers (Ci.i = C2,i = ... = Cm,i) and may be computed only once. The 
different SF's are applied to the Dedicated Physical Data Channel (DPDCH) and the 
M Rake receivers operate on the data channel. In 3GPP a single DPDCH is spread 
by the code Cd,i =Cch.sF,k where k=SF/4. A consequence of this specification is that 
the code employed in the DPDCH is known once the SF is known. 

Despreading with a hypothetical SFi that is smaller than the actual SF2 
produces a decision statistic that is part of the actual one. Accumulating the decision 
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statistic for SFi over SF2/SF1 intervals produces the actual statistic obtained by 
despreading with the code corresponding to SF2. This property allows despreading 
to occur with SF=4 for all possible SF's and the partial despreading results to be 
combined once the actual SF is known. 

Despreading with a hypothetical SFi which is larger than the actual SF2 will 
produce a decision statistic that nnay contain all or parts of the SF1/SF2 decision 
statistics for the actual SF2 depending on the combination of SF1/SF2 data symbols 
for SF2. If all SF1/SF2 data symbols for SF2 are identical, then the decision statistic 
obtained after despreading with SFi will be identical with the sum of SF1/SF2 decision 
statistics obtained by despreading with SF2. If those data symbols are not identical, 
cancellations will occur during the despreading with SFi and the absolute value of the 
resulting decision statistic will be smaller than the one obtained by summing the 
SF1/SF2 absolute values obtained by despreading with SF2. This property of 
statistical averaging will be exploited by the preferred embodiment SF estimation 
methods. Its validity relies on the fact that the decision statistics (symbols from MRC 
for the SF hypotheses) will be obtained over several slots. 

Based on these observations, preferred embodiment SF estimation methods 
are based on the statistical properties of the resulting decision statistics. SF 
estimation thus becomes a detection problem with M hypotheses (M < 7). For an 
additive white Gaussian noise (AWGN) channel, SF estimation for minimum error 
probability is based on the likelihood ratio test (LRT). For normalized variance <t^, the 
LRT simply chooses the largest of the likelihood metrics 

Li =ln(Pi) - (ISil - ^i)2/(2cy2), i=1 ,2 M, (1) 

where ln(.) is the natural logarithm, Pi is the a priori probability for hypothesis i (i.e., 
the hypothesis that the actual spreading factor equals SF|), S| is the output from the 
i*'' Rake receiver using MRC, m is the mean, and is the common variance (after 
some normalization). If the a priori probabilities are the same for all possible SF's, 
the ln(P|) terms can be omitted from the likelihood metrics. Also, because the 
DPDCH data channel is modulated, use the absolute value |Si| of the Rake output. 
Although, by taking the absolute value of the Rake output, the resulting distribution is 
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no longer Gaussian (it is Rician), the LRT in (1) is still applicable to the SF estimation 
problem as it will become apparent In the following. 

Because the typical signal-to-nolse ratio (SNR) of Eg/No in 3GPP ranges 
between 0 and 6 dB, after one symbol for the correct hypothesis (i.e., SF value) 
^i^/2cy^ ranges between 0.5 and 2 due to im^=EJ2 and 2a^=No. The probability Pf of 
incon-ect SF estimation for such SNR values is much larger than the typical FER 
values in 3GPP. Moreover, because of fading and power control errors, even for 
much larger SNR values, Pf can be prohibitively large if the decision statistic is 
obtained after only one despread symbol. Therefore, the absolute values of the 
Rake outputs should be accumulated over a number, P, of symbol periods. 

If Si is a Gaussian-distributed random variable with mean \i\ and variance a-^, 
then |Sj| has a Rice distribution with mean 

E[ISil] = i^i erf(|ii N2 ai) + V(2/7i) Cj exp(-ii{^ 12 o^) 
and variance 

Var[ISil] = Oi2 + ^i2-E[ISil]2. 
The error function, erf(.), is defined as erf(x) = {2Hn) f[o.x] exp(-y2)dy. Notice that 
when \i?l2cy? > 2 (i.e., SNR > 3dB), then E[ISil] ^ ^ij and Var[ISil] ^ ai^. However, 
in 3GPP, because of fading and power control errors, the SNR may fall considerably 
below 3 dB and the previous approximation will not be valid in general, and 
especially for low SNR's where blind SF estimation is most difficult. 

Decision statistics (DS) will be obtained by summing the absolute values of the 
Rake outputs, ISjl, over P symbol periods. Thus, the decision statistic for the ith 
hypothesis, DS, = (1/P) ElSil. The central limit theorem (CLT) implies each DS\ is 
Gaussian distributed. The CLT is applicable because P will typically exceed 40 
(smallest number of symbols per slot is 10 at the value for SF=256, and the 
observation period is 4 slots). Its distribution (pdf) will therefore be completely known 
if its mean and variance are known. The next step then is to find the mean and 
variance of the decision statistics for each of the M hypotheses. As was previously 
noted, the despreading result for each hypothesized SF depends on whether the 
actual SF is larger or smaller that the hypothesized SF. Therefore, the mean and 
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variance of the decision statistics will depend on the actual SF as detailed in the 
following sections. 

3. Decision Statistic Mean and Variance for actual SF larger than hypothesized SF 

First consider the case of when the actual SF2 is larger than the hypothesized 
SFi. While accunfiulating the DS2 over P SF2 symbol periods, DSi is accumulated 
over (SF2/SFi)P SFi symbol periods. Since H2= (SF2/SFi)|j,i and ctsf2 = 
(SF2/SFi)^'Wi. the mean of DS2 is 

E[IS2|] = }i2 erf(fi2 /V2 CT2) + V(2/7r) 02 exp(-^2^ /2 02^) 
Normalizing with V2 02 yields: 

E[IS2|] = VSNR erf(VSNR) + VI/ti exp(-SNR) 
and the mean of DSi is 

E[ISi|] = {SFi/SF2)n2 erf(^2 V(SFi/SF2)/V2 C2) 

+ V(SFi/SF2)V(2/7t)cy2 exp(-H2^ (SFi/SF2)/2 az^) 
Multiplying E[ISil] by SF2/SF1 to account for the longer averaging for DSi yields 
E[ISi|] = [12 erf(|Li2 V(SFi/SF2)/V2 ctz) 

+ V(SF2/SFi)V(2/7t) G2 exp(-p.2^ (SFi/SF2)/2 (J2^) 

or 

E[ISi|] = VSNR erf(V(SFi/SF2)VSNR) + V(SFi/SF2)Vl/7t exp(-(SFi/SF2)SNR) 
In general, the difference between the means increases as the SNR increases and 
as the ratio of the spreading factors increases. 
The variance of DS2 is 

Var[IS2l] =cy2^ + ^2^ - E[IS2l]^ 

= cs2^ + )x.2 - [pi2 erf(n2 /V2 csi) + V(2/7t) ©2 exp(-^2^ /2 <52)f 
The variance of DSi is 

Var[ISil] = (SFi/SF2)a22 + (SFi/SF2)V2^ - E[ISilf . 
Since DSi is accumulated SF2/SF1 times more than DS2 and the result is scaled by 
(SFi/SF2)^^^ the variance of DSi is 

Var[ISil] = 02^ + - [^2 erf(V(SFi/SF2)Ki2/V2 G2) 

+ V(SFi/SF2)V(2/7t)a2 exp(-(SFi/SF2)^i2^/2CT2^)]2 
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It should be noted that since SF2 is not known, the normalizing factor V2 02 is 
also not known. For that reason, a reference variance computed over a 256-chip 
interval can be used. The effect of this is simply to progressively divide the decision 
statistics by powers of V2. Thus, the DS for SF=256 is divided by 1.0, the DS for 
SF=128 is divided by V2, the DS for SF=64 is divided by 2, the DS for SF=32 is 
divided by 2V2, the DS for SF=16 is divided by 4, the DS for SF=8 is divided by 4V2, 
and the DS for SF=4 is divided by 8. In general, the decision statistics are in effect 

divided by (V2)'-^ where SF=2', 1=2,3 8. Notice that this nomializatlon takes place 

automatically when dividing with V2 <J256- The reason for dividing by V2 g2 is to make 
the decision statistics dependent on the SNR and thus to be able to know their 
theoretical values based on the measured SNR. Measuring the variance 02 is 
needed for power control purposes in 3GPP and it is not a computation introduced by 
the preferred embodiment SF estimation. Also, as shown in the following, this 
normalization is not necessary for the preferred embodiment SF estimation methods. 

4, Decision Statistic Mean and Variance for actual SF smaller than hypothesized SF 

In contrast to the foregoing section, now consider the case of actual SF2 
smaller than hypothesized SFi . The value of the decision statistic for the 
hypothesized SFi depends on the transmitted data symbols since despreading with 
SFi occurs over several SF2 symbols. Since the absolute value of the despread 
signal is taken, the mean of DSi does not depend on the sign of the data symbols, 
but only on the number of positive and negative ones. Denote SF1/SF2 as N, then 
the probability that / out of N transmitted SF2 symbols will be different is expressed in 
terms with the binomial coefficient as 

2 r 72'^ 

where the extra factor of 2 reflects the fact that negative and positive are equivalent. 
The corresponding mean is iN-2i)n2. Therefore, 

E[|Si|] = 1/2'^-^ [ 
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{iV^2.erf(A/ uzHl^Ncjz) + VA/V(2/7t) 02 exp{-l^\i2^l2Na2^)} 



{(A/-2)p2.erf((A/-2)^i2/A/2VA/a2) + VA/V(2/7i) Q2 exp{-{N-2f |X2^/2A/<J2^)} 



f N ^ 
,iV/2-l 



{2)X2.erf(2p2/V2VA/cT2) + VA/V(2/7i) 02 exp(-2 \i2^l2Nc5^)} 
{1/2 V/VV(2/7t) a2} ] 



or by normalizing witli V2a2 
E[|Si|] = Ml"-' [ 



N 



0 ; 



{A/VSNR erf(A/VSNR) + <N-i{Mn) exp{-N SNR)} 



.1/ 



{(/V-2)VSNR erf((A/-2)VSNR /V/V) + VA/V(1/7c) a2exp(-(A/-2)2SNR/VA0} 



f N 



' N 



Nil, 



{2VSNR erf(2VSNR HN) + VA/V(1/7i) exp(-22 SNR/A/)} 
{1/2 VA/V(1/jt) } ] 



The mean for DS2 after accumulation over N symbol periods is equal to N times the 
mean over one symbol period, and because of the scaling by (1/A/)^'^ (for noise mean 
component normalization), the mean is 

VA/ E[IS2|] = VA/ h2 erf(p.2 /V2 02) + VA/V(2/7t) 02 exp(-n2^ 12 
or, by normalizing with V2 02: 
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Va/ E[IS2|] = VA/VSNR erf(VSNR) + VA/V(1/7r) exp(-SNR) 
The variance of DS2, after accumulating over N symbol periods and normalizing, is 

Var[IS2l] = 02^ + - E[\S2\f. 
The variance of DSi is 



Var[ISil] = ^2^/2'^"^ S 



(A/-2/)2 + A/cj2^-E[ISilf. 



5. Estimation methods 

The larger variance for the incorrect hypotheses will not be an issue for SF 
estimation because DS is ultimately obtained over hundreds of symbols (for an 
observation period of 4 slots, the smallest number of symbols is 10 per slot when 
SF=256 and thus totals 40 symbols). It is the mean value of the various hypotheses 
for each of the M=7 possible SF values that will determine the detection strategy. 
Since the mean DS value of each hypothesis depends on the correct SF and the 
hypothetical SF being tested, the mean values need to be computed for each SF 
separately. The post-processing of the averaged absolute values after MRC for each 
possible SF should be such that the final DS for the correct SF hypothesis achieves 
maximum distance from the final OS's for all incorrect SF hypotheses and the 
distance among incorrect hypotheses is practically the same. 

The objective of any SF estimation method is to ensure that the probability of 
incorrect SF estimation is substantially smaller than the corresponding FER and 
preferably smaller than the TFCI codeword decoding error rate. Also, SF estimation 
should be perfomied in the shortest possible observation period in order to maximize 
its benefits. Clearly, the observation period cannot exceed 15 slots, or 1 frame 
period, because the SF can then be determined from the TFCI. The maximum 
observation period determines the number of samples that need to be stored before 
despreading and the goal is to limit this value to 50% or less of the frame period. The 
observation period depends on the environment, due to the dependence of the DS 
reliability on the SNR variations, and on the actual SF's to be detemiined because 
this specifies the number of symbols in the observation period. Large SF's (e.g. 256, 
128) require long observation periods because of the corresponding low data rate. 
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Small SFs (e.g. 4, 8, 16) require shorter observation periods because of the 
corresponding high data rate. If the minimization of the observation period is 
preferable to the estimation of each possible SF value, the SF estimation process 
can be used to determine the maximum SF value that can be reliably decided for the 
specified observation period. For example, assume that for an observation period of 
4 slots the SF detennination can be performed with the required reliability for 
SF=4,8,16,32, and 64 but not for SF=128 and 256. Then, if the SF is determined to 
be 64 or larger, despreading may proceed with SF=64. If at the end of the frame the 
actual SF is determined to be 128 or 256, the MRC results based on SF=64 can be 
appropriately combined. This implies that the MRC results need to be stored if the SF 
is determined to be larger than 64, but the size of the memory is smaller by a factor 
of 16 relative to the always possible option of despreading with SF=4 and storing the 
con-esponding MRC results for combining at the end of the frame. 

There are several methods for post-processing of the decision statistics, each 
defining a corresponding preferred embodiment SF estimation method. The 
objective is to ensure the equivalence of incorrect decision statistics while achieving 
a larger metric for the correct one. Figures 2-3 illustrate the common steps for 
preferred embodiment SF estimation methods and include the steps of: (a) 
Despreading with multiple SF's and performing MRC; (b) Taking the absolute value 
of MRC output; (c) Weighting the absolute value with the SNR estimate at the 
corresponding time period; and (d) Accumulating the weighted decision statistics over 
the entire observation period. Notice that because the Walsh code for each possible 
spreading factor is simply a repetition of the Walsh code for SF=4, all statistics can 
be obtained with only one despreading process. Two preferred embodiment 
post-processing methods for the decision statistics obtained after step (d) will be 
described in detail. The first is based on the ratios of the decision statistics while the 
second is based on a linear combination of them. 

6. First preferred embodiment 

The first preferred embodiment method evaluates the ratio R(SFi) = 
D(SFi)/D(SFi/2) where D(SF) is the decision statistic DS after step (d) for SF. The 
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ratios R(SFi) for M=7 hypothesized SF's for each actual SF are shown in Figures 
4,1-4.7. In particular, Figure 4.1 shows the case of the actual SF = 256 and the 
ratios R(SFi) with SFi = 4, 8, 16, 32, 64, 128, 256, and 512 where estimates for 
D(512) and D(2) were used in R(512) and R(4), respectively, instead of actual 
despreading with SF=512 and SF=2. Similarly, Figure 4.2 shows the same ratios but 
for the case of the actual SF = 128, and so forth down to Figure 4.7 which shows the 
same ratios for the case of the actual SF = 4. 

Several observations can be made from Figures 4.1-4.7. First, for SFi equal 
the actual SF, the ratio R(SFi) is the largest and R(2SFi) Is the smallest. Second, for 
SF's smaller than the actual SFi, the value of R(SF) decreases as the SF decreases 
while for SF's larger than the actual SFi, the value of R(SF) remains approximately 
the same and equal to 1/V2 for the SNR's in the range of interest. This reflects the 
statistical averaging of the data symbols for SF's larger than the actual one. Third, 
the ratio R(SF) for actual SFi equals R(2SF) for actual 2SFi. Fourth, the difference 
R(SFi)-R(2SFi) is positive and the largest for SFi the actual SF, while for other SF's, 
the difference R(SF)-R(2SF) is either negative or zero. 

The decisions for the first preferred embodiment method using R(SF) are 
based on the previous observations. In addition to R(SF), the differences 
R(SF)-R(2SF) are also computed and serve as the primary post-processing statistic 
for other first preferred embodiment methods. The nature of R(SF)-R(2SF) is similar 
to the statistic for orthogonal signaling. This is most clearly observed for SFi=8 
where all other differences are essentially equal to zero while the difference 
R(8)-R(16) is positive. If R(SF)-R(2SF) is not zero, it is negative (a small negative 
value compared with the positive value for the correct SFi) and this in effect adds to 
the R(SFi)-R(2SFi) for the actual SFi. The differences of the ratios for the various 
SF's are presented in Figures 5.1-5.7. 

In order to have a complete set for the differences R(SF)-R(2SF), the statistics 
for SF=512 and SF=2 need also be obtained. However, simplifications can apply for 
both SF=512 and SF=2 as follows. 

For all SF's and the SNR range of interest (<10 dB), R(512) is practically equal 
to 1/V2 except for actual SFi=256. Even then, R(512) becomes considerably smaller 
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than I/V2 only for large SNR values where the decision statistics are reliable and the 
SF estimation can be made with R(256)-R(512) being smaller by 0.0-0.05 compared 
to its actual value. A fixed value for R(512) is also beneficial because the 
computation of R(512) in the case of actual SFi=256 would have low SNR and thus 
large variance. Figures 4.1-4.7 and 5.1-5.7 assume a value of R(512) equal to 1/V2 . 

A simplification similar to the one for SF=512 can be made for SF=2. Based 
on the SNR measurements during a particular slot, the value of R(2) can be set to its 
theoretical one for SF=4. The decisions for the SF will be based on the difference 
between the R(SF)'s and the effect of the approximation for SF=2 is to benefit the 
decision for every SF other than SF=4. For SF=4 the approximation provides results 
that are practically equivalent to the actual ones. Figures 4.1-4.7 and 5.1-5.7 assume 
this approximation for R(2). 

The first preferred embodiment methods' decision for the SF can be based on 
only the difference between the R(SF)'s as previously described. Alternatively, more 
elaborate schemes can be devised based on both the differences between the 
R(SF)'s of Figures 5.1-5.7 and the actual R(SF) values of Figures 4.1-4.7. 

More particularly, a first preferred embodiment blind SF estimation system for 
3GPP uplink in FDD mode could proceed as follows. (1) At the start of a frame in the 
data channel (DPDCH) received from a user, begin applying the scrambling code and 
the spreading code (repeating +1 +1 -1 -1) in a Rake detector, and at intervals of 2' 
chips take the accumulated sum as an Sj value, then take the magnitudes |Sj| of the 
soft outputs to start generating the various decision statistics DS(2'). (2) Once P 
symbols for SF = 2' have been received (that is, after P intervals of 2' chips in which 
the intervals may be separated by unused intervals), compute the ratio R(2') = 
DS(2')/DS(2'"^) and compare this to the other ratios already computed (which 
correspond to smaller powers of 2) for a decision opportunity: if R(2*"^) - R(2') 
exceeds a positive threshold (such as 0.1), then decide the spreading factor for the 
frame as SF = 2*"^ whereas if R(2'"^) - R(2') is about 0 or negative, then the SF 
decision is deferred until R(2'*^) is tested; and so forth. See Figures 5.1-5.7. Note 
that when R(2'"^) - R(2') is about 0 or negative and more accumulating is occurring, 
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the previously computed DS(4), DS(8), DS(2') may be recomputed to increase 
accuracy of the ratio differences R(2'"^) - R(2'). 

Alternatively, in addition to R(2'"^) - R(2') exceeding a positive threshold, a 
decision can include testing R(2') to be at most about 1/V2. 

Further first preferred embodiments proceed as follows. Initially, find two 
candidate SF's, SF2 and SF3, by taking SF2 so that R(SF2) - R(2SF2) is the largest 
ratio difference and SF3 so that R(SF3) - R(2SF3) is the second largest. Next, decide 
SF to be SF2 if R(SF2) is at least a threshold percentage (e.g., 7-10%) larger than 
R(SF3) and conversely to be SF3 if R(SF3) is at least a second threshold percentage 
(e.g., 7-10%) larger than R(SF2). However, if the ratios are within threshold 
percentages of each other, then decide SF to be SF2 when R(SF2) - R(2SF2) is at 
least another threshold percentage (e.g, 8-12%) larger than R(SF3) - R(2SF3). Then 
if none of the previous decision conditions holds, and if either R(SF2) or R(SF3) is a 
further threshold percentage (e.g., 2-3%) larger than the other, then take SF to the 
one with the larger ratio. Lastly, if no decision has been made by the foregoing, take 
SF = SF2. 

7. Second preferred embodiments 

The second prefenred embodiment methods post-process the decision 
statistics according to the formula 

F(SF)=F(i)=[DS(i) - aDS(i+1)] / p^"' 
where SF=2', i=2,3,...,8 and a, p are real numbers. Using similar arguments as with 
the observations in the first preferred embodiment methods, the factor a is selected 
equal to V2. The result of this will be that if the actual SF is smaller than the ones 
considered in the above expression, F(SF) will be zero. The factor p can also be 
viewed as a weighting factor, trying to compensate among the different levels of 
reliability for the decision statistics. A value of p larger than 1 .0 assists the decision 
statistics for the larger SF's while a value of p smaller than 1 .0 assists the decision 
statistics for the smaller SF's. 
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As with the first preferred embodiment methods, the values for SF=512 and 
SF=2 are also needed. Using similar arguments, despreading for SF=512 can be 
avoided and F(256) can be set to 0.0. For SF=2, despreading can also be avoided 
by setting SF(2) to its theoretical value, for the measured SNR during the particular 
slot, for a hypothetical SF=4. This will produce accurate results when the actual 
SF=4. For all other values of the actual SF, F(2)=DS(2)-DS(4) will be larger than its 
actual value and will typically be the largest value. This is, however, beneficial for the 
final decisions which will be based on F(SF)-F(SF/2). Thus, the previous setting for 
F(2) will actually improve the estimation of all SF's greater than 4, and will produce 
accurate results for SF=4. Figures 6.1-6.7 and 7.1-7.7 better illustrate this point. 
Figures 6.1-6.7 show the decision statistic for a =V2 and p = 1 . 1 . The curves for F(2) 
and F(256) were based on the previously discussed approximation. The differences 
between decision statistics and their properties resemble the ones for the first 
method and similar observations apply. Notice that except for SF=4 and SF=8, F(2) is 
the largest for all other SF's with the second largest value attained for the actual SF. 
As with the first method, the decision for the SF can be based only on the difference 
between the F(SF)'s. Alternatively, both the difference between F(SF)'s (Figures 
7.1-7.7) and the F(SF)'s may be considered. Because of the much larger reliability of 
the decision statistics for the smaller SF's as a result of more data symbols, the 
decision statistics for the smaller SF's do not need to be obtained for every symbol 
during the observation period. For example, for an observation period of 8 slots, the 
decision statistics for SF=256, SF=128, and SF=64 are obtained during every slot, for 
SF=32 and SF=16 they may be obtained only every other slot, while SF=8, and SF=4 
they may be obtained during only 2 slots. This considerably simplifies the 
computational complexity of the SF estimation method since MRC for the smaller 
SF's need only be perfonned 1/4 or 1/2 of the time. 

8. Altemative preferred embodiments 

Both first and second preferred embodiment SF estimation methods may be 
used simultaneously. Because of the different statistical properties of the 
con-esponding decision statistics, the two methods may produce different SF 
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decisions. If despreading occurs for both decided SF's, this can be used to reduce 
the overall SF estimation error rate. The final SF decision will be based on the TFCI 
and the incorrect despreading results will be dismissed. This implies that the MRC 
results for the smaller of the two different SFs are stored. 

Because of the lower reliability and larger variance of the decision statistics 
during a given observation period, for SF=256 and SF=128, the resulting decisions 
are the least reliable and despreading for both SF=256 and SF=128 may be 
performed if either is decided. Using similar reasoning, if SF estimation errors for any 
particular SF are mostly caused by another specific SF (e.g. actual SF of 256 and SF 
causing most errors of 128), then despreading may occur for both such SF's. Again, 
in such cases, the final SF decision will be based on the TFCI assuming that the 
MRC results for the smaller of the two SFs are stored 

Likewise, the high data rate situation (multiple data channels) with more than 
one SF = 4 code can be handled by despreading with all 3 of the possible SF = 4 
codes (+1+1-1-1, +1-1+1-1, and +1-1-1+1) when the decision is SF = 4, and use the 
TFCI at the end to determine which of the results to use and which to discard. 
Another alternative is to despread with each of the possible SF = 4 codes over a few 
slots, measure the SNR, and determine the code existence in that manner. Also 
notice that in the 3GPP uplink, compressed mode is supported only through SF 
reduction by a factor of 2 (no puncturing) and this poses no problem for SF 
estimation. 

9. Experimental results 

Simulations with statistics accumulated over 8 slots and over 6 slots (i.e., 
varying P), with mobile speeds of 3, 30, and 100 km/h, Rayleigh fading, 1-2 
multipaths, and 5,7,or 10 uncoded bits generating the 30 TFCI coded bits yielded the 
following. 

For an observation period of 8 slots, the preferred embodiment SF estimation 
methods have an enror rate that is substantially smaller than the 10 or 7 uncoded bit 
TFCI decoding enror rate even for the SF's of 256 and 128. The SF estimation error 
rate is in general comparable with the decoding error rate for a TFCI with 5 uncoded 
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bits. Also, an observation period of 6 slots or smaller will generally suffice for the 
smaller SF's, and it will suffice for all SF's in slow Rayleigh fading, rich multipath 
environments, base stations with antenna diversity, or AWGN channels. 

The SF estimation error rate for SF = 256 and SF = 128 will significantly 
decrease if despreading for either SF = 128 or SF = 256 continues after a respective 
decision for SF = 256 or SF = 128 is made. This is because the vast majority of 
errors for SF = 256 are caused by deciding SF = 128, while for SF = 128 by deciding 
SF = 256. By despreading both SF's a final decision can be made at a later instant 
or at the end of the frame. The penalty for despreading with both SF = 128 and SF = 
256 is the additional memory required to store the extra symbols. This memory and 
the additional processing (MRC) are small because of the low data rate for these two 
SF's. An additional benefit is that the observation period may then be decreased to 
about 6 slots or less. A compromise between longer observation periods and SF 
estimation may be to determine whether the SF is smaller than 64. If it is, the actual 
SF is determined. If it is not, despreading can be performed with SF = 64 and the 
results combined once the actual SF becomes known at the end of the frame. This is 
the same as despreading with SF = 64 and combining the results at the end of the 
frame, only now the required processing and memory are much smaller. 

In short, the preferred embodiment blind SF estimation methods are viable in 
3GPP because their reliability is better than the frame error rate and comparable to 
the TFCI decoding error rate. The associated complexity is small because the 
operations needed are a single despreading, MRC, computation of the absolute MRC 
values, and accumulation of the absolute values over several slots. A DSP can 
process the final values and make an SF decision. An observation period of 5-8 slots 
produces the desired results. This translates to a reduction of about 50-67% in the 
buffer size and data transfer rate as compared to relying upon TFCI decoding SF 
determination. The observation period can be reduced to 5 slots or less if both SF = 
128 and SF = 256 are despread when a decision for either is make. Further 
reductions may also be possible for base stations employing antenna diversity. 
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10. Modifications 

Several modifications of the preferred embodiment SF estimation metfiods 
may be devised while retaining the feature of steps like (a) Despreading with all 
possible SF's (and performing MRC of resolved multipaths); (b) Taking the absolute 
values of (MRC) outputs; (c) Weighting the absolute values with the SNR estimate for 
the corresponding time period; (d) Accumulating the weighted absolute values to 
generate decision statistics; and (e) Use the decision statistics for the various 
spreading factors to decide on an SF estimate. Step (c) is also optional. 

For example, various quantities can be changed, such as P (the number of 
symbols averaged for DS); M (the number of possible spreading factors); the number 
of slots over which the (MRC) outputs are accumulated which may adapt to the SNR; 
other formats with differing number of chips per slot, number of slots per frame; 
different spreading code structures with the case of the codes for smaller SF's 
imbedded in the codes of larger SF's permitting the use of a single Rake detector; 
and so forth. 
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